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Barlowite Cui{OH)c,FBr is a newly found mineral containing Cu^^ kagome planes. Despite 
similarities in many aspects to Herbertsmithite Cu 3 Zn{OH)gCl 2 , the well-known quantnm spin 
liqnid (QSL) candidate, intrinsic Barlowite turns out not to be a QSL, possibly due to the presence 
of Cu^+ ions in between kagome planes that induce interkagome magnetic interaction [PRL, 113, 
227203 (2014)]. Using first-principles calculation, we systematically study the feasibility of selective 
substitution of the interkagome Cu ions with isovalent nonmagnetic ions. Unlike previous speculation 
of using larger dopants, such as Cd?'^ and Co?'^, we identify the most ideal stoichiometric doping 
elements to be Mg and Zn in forming CuiMg(OH)^FBr and Cu3Zn{OH)eFBr with the highest 
site selectivity and smallest lattice distortion. The equilibirium anti-site disorder in Mg/Zn- doped 
Barlowite is estimated to be one order of magnitude lower than that in Herbertsmithite. The single¬ 
electron band structure and orbital component analysis show that the proposed selective doping 
effectively mitigates the difference between Barlowite and Herbertsmithite. 


Quantum spin liquid (QSL) represents a new state of 
matter characterized by long-range entanglement, be¬ 
yond the conventional symmetry-breaking paradigirf^. 
Realizing QSL in real-world materials has been a long- 
sought goal for decaded^m The most promising candi¬ 
date so far is Herbertsmithite Cu 3 Zn{OH)eCl 2 , which 
realizes the S' = 1/2 antiferromagnetic (AFM) Heisen¬ 
berg model on the 2D kagome lattic^. Extensive the¬ 
oretical studies have suggested that this model is likely 
to achieve a QSL ground-state, despite close in energy 
with other competing phases^n^. Experiments on Her¬ 
bertsmithite have also shown QSL-like features, such as 
the abse nce of any observed magnetic order down to 50 
]^pj i4 | i5 | unusual continuum of spin excitationJi^. 

However, the inevitable Cu/Zn antisite disorder makes 
the interpretation of experimental data difficulffi^. It re¬ 
mains an open debate whether these defects obscure the 
intrinsic signals, such as a tiny spin gap that is crucial 
for the classification of the ground state^S. 

Very recently, Barlowite Cui(OH)c,FBr as a new 
kagome compound was discoverecpSI. its structure closely 
resembles Herbertsmithite, whereas the Cu/Zn antisite 
disorder is automatically avoided. Therefore, studies on 
this new material are expected to shed fresh light on un¬ 
derstanding the kagome physics and QSL phase. Inter¬ 
estingly, Barlowite is diagnosed with a Curie-Weiss con¬ 
stant 0 c' 14/=-I36K close to Herbertsmithite, yet it un¬ 
dergoes a spin-ordering phase transition at I5KP^. The 
low-temperature magnetic properties were further in¬ 
vestigated by combining first-principles calculation with 
experiments^il. Since the main structural difference be¬ 
tween these two materials is the cations occupying the 
interkagome sites, i.e. and in Barlowite and 

Herbertsmithite, respectively, it is suggested that sub¬ 
stituting the interkagome sites with nonmagnetic ions 


should tune Barlowite into the same phase as Herbert¬ 
smithite. Specifically, relatively larger elements, such as 
Sn and Cd, are speculated as possible candidates for sub¬ 
stitution based on the simple argument of lattice spacing 
of the interkagome sites^. 

Cu2->M2+(S=0) 



FIG. 1: Unit cell and atomic structure of Barlowite, and the 
schematic illustration of the proposed selective doping. 

In this Letter, we identify the most promising candi¬ 
dates for realizing the selective doping to form stoichio¬ 
metric doped Barlowite^based on density functional the¬ 
ory (DFT) calculation^. We systematically calculate 
the doping energies and analyze the doping selectivity of 
a series of nonmagnetic group 2 and 12 elements. Unlike 
the previous speculatioiP^, larger dopants are found to 
have lower site selectivity and tend to distort the kagome 
plane more than smaller dopants. Most importantly, 
we identify Mg and Zn to be the most ideal choices of 
dopants to form stoichiometric Cu^Mg{OH)QFBr and 
Cu 3 Zn(OH)c,FBr compounds in the Barlowite family, 
with the highest site selectivity in substituting the interk- 
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agome Cu ions and the least lattice distortion in kagome 
planes. Statistical analysis shows that the equilibrium 
distribution of Mg/Zn in the Mg/Zn-doped Barlowite at 
the typical growth temperature exhibits a level of an¬ 
tisite disorder significantly lower than that in Herbert- 
smithite. Also, single-electron band structures of intrin¬ 
sic and doped Barlowite are calculated, and discussed in 
comparison with Herbertsmithite. 

Figure 1 shows the atomic structure of Barlowite. Sim¬ 
ilar to Herbertsmithite, it contains Cv?'^ kagome planes 
connected by hydroxyls. The difference lies in the in- 
terkagome site: in Barlowite, there are additional Cu^'^ 
ions between the kagome planes, which act as additional 
spin 1/2 centers and mediate interkagome spin exchange. 
Therefore, to clarify the different magnetic ground states 
between Barlowite and Herbertsmithite, one way is to 
remove these out-of-plane spins by selective doping. The 
chosen dopants should be spin zero and isovalent. Using 
these two criteria, we have considered elements of Mg^ 
Ca, Sr, Ba from group 2 and Zn, Cd from group 12. The 
ionic radius (ri)and standard chemical potential in aque¬ 
ous solution(^/) of these ions are listed in Tab. I. Note 
that and Zn^~^ are close in radius to Cu^^, while 

the other ions are larger. In addition, all these dopants 
under consideration have a smaller than CvS''^, indi¬ 
cating a stronger tendency to stay in the solution. 


TABLE I: Ionic radius (vi) and standard chemical potential 
in aqueous solution (^/) of Cu and the dopants under consid¬ 
eration 



n (pmj^ 

(eVj^ 


73 

0.68 


72 

-4.73 

Ca^+ 

100 

-5.73 

Sr 2+ 

118 

-5.78 


135 

-5.81 


74 

-1.52 

Cd2+ 

95 

-0.80 


The doping process is expected to take place by adding 
dopant ions in solution during the hydrothermal growth. 
The net reaction equation is written as : 

{l-Xd)Cui{OH)QFBr + AxdM'^+ + QxdOH- 
+XdF~ +XdBr~ [Cui-a:^M^^]4{OH)eFBr 


in which denotes the dopants. 

Equation (1) can be considered as a combination of 
two subreactions: 


Cu4{OH)eFBr + 4:XdM^+ ^ 
[Cui.,,,M,MOH)eFBr + AxdCu'^^ 


( 2 ) 


AxdCu^'^' + QxdOH + XdF + XdBr 
—>■ XdCu 4 {OF[)QF Br 

Equation (2) describes the simple substitution process. 
Then, after the Cu ion is exchanged into the solution. 


the overall free energy of the system can be lowered by 
forming more deposits of Barlowite [Eq.(3)]. The latter 
provides the thermodynamical driving force for the dop¬ 
ing process to proceed. Equation (3) is actually nothing 
but the growth of undo ped Barlowite as reported in pre¬ 
vious experiment Hence, we will focus on evaluat¬ 

ing the experimental feasibility of Eq.(2) only. 

The central physical quantity we are going to calcu¬ 
late is the standard doping energy Ed as defined by the 
total energy difference per substitution. According to 
Eq.(2), Ed consists of two parts: Ed = AAb-I-A/x/, where 

AAb = {E[Cui-,^^M,,Ji{OH)eFBr — Ecui,(OH)fiFBr)I^Xd is 

the energy change of Barlowite after doping and A/i/ = 
/ip^ 2 + — standard chemical potential differ¬ 

ence between the two ions (Tab. I). There are two in¬ 
equivalent doping sites: Cul is in the kagome plane; Cu2 
is between the kagome planes. We use Edi and Ed 2 to 
differentiate these two types of doping energies. Their 
difference AEd{xd) = Edi{xd) — Ed 2 {xd) tells the site 
preference for dopants, i.e. defines the degree of selective 
doping. 
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FIG. 2: (a) Ed and (b) AEd of different dopants. E^ in (a) 

corresponds to the standard doping energy for the growth of 
Herbertsmithite. The dotted line in (b) only serves as a guide 
to eyes, (c) Equilibrium anti-site disorder of Herbertsmithite 
(H) and Zn/Mg doped Barlowite (B-Zn/Mg) 







































3 


Our calculation on /S.Eb is carried out using the 
VASP packagd^, which solves the DFT Hamiltonian self- 
consistently using the plane wave basis together with the 
projector augmented wave methocP^. A plane-wave cut¬ 
off of 500 eV is enforced. The self-consistent iterations 
are converged to 0.1 meV precision of the total energy. 
We employ the standard generalized gradient approxi¬ 
mation (GGA) for the exchange-correlation functional, 
which is known to satisfactorily describe ionic bond¬ 
ing and cohesive energy even for transition element^^. 
Barlow ite st ays in the paramagnetic phase above 15K 
(lmeV)2212iJ. The Dzyaloshinskii-Moriya (DM) interac¬ 
tion arising from spin-orbit coupling is estimated to be 
of the order of 1 me\^. Therefore, a spineless calcula¬ 
tion on the structure and total energy is reasonable. As a 
benchmark. Table 2 summarizes the calculated structural 
parameters of undoped Barlowite, which agree with the 
experimental results well. Note that we do not intend to 
discuss strong correlation effects associated with the Cu 
3d orbitals within this methodology. 

To simulate doping, we construct a 2 x 2 x 1 supercell 
containing 24 in-plane Cul sites and 8 interkagome Cu2 
sites, and replace one of the Cu with a dopant. We first 
fix lattice contants to experimental values, while let the 
atomic coordinations fully relax until the forces are less 
than 0.05 eV/A. The total-energy integration over the 
Brillouin zone is obtained on a T-centered 2 x 2 x 2 k- 
mesh. 


TABLE B: A comparison of structural parameters between 
theory and experiment 


Barlowite 

Exp.l^ 

Exp.l^ 

Cal. 

Lattice a/b (A) 

6.68 

6.80 

6.73 

Lattice c (A) 

9.31 

9.31 

9.47 

Angle Cul-O-Cul 

117.4° 

117° 

117.3° 

Angle Cul-0-Cu2 

95.8° 


96.6° 


Figure 2(a) shows the calculated Ed of different 
dopants occupying the Cul and Cu2 sites. As 
a reference, we have also calculated the standard 
doping energy for the growth of Herbertsmithite 
[E^ in Fig. 2(a)] as described by the follow¬ 
ing equation: [Cuz+xaZni_xj\i.OH)QCl 2 + XdZn^+ ^ 
Cu'iZn{OE[)^Cl 2 + XdCu^^. The positive doping en¬ 
ergy acts as a reaction barrier, which limits the kinet¬ 
ics of Eq.(l). The value of Ed 2 for Zn is found almost 
the same as again reflecting the similarities between 
Barlowite and Herbertsmithite. Ed typically increases 
with the ionic radius as a size effect: Mg is even easier 
to substitute Cu than Zn, while larger dopants are more 
difficult. 

The kagome-site doping energy {Edi) is always higher 
than the inter kagome-site doping energy {Ed 2 ), indicat¬ 
ing that the latter is the preferred site for doping. The 
preference for dopants to occupy the interkagome site 
provides exactly the type of doping selectivity we need. 
It can be understood by noticing the larger space around 


the Cu2 site than the Cul site (Fig.l). In Fig. 2(b), we 
plot AEd = Edi — Ed 2 as a function of ion radius, show¬ 
ing directly the site preference in relation to ionic radius. 
Within each group, AEd decreases monotonously. The 
elements from two groups are however not on the same 
line. This result is reasonable considering that the out¬ 
ermost shells of group 2 and 12 ions are p-electrons and 
d-electrons, respectively. Thus, they will exhibit differ¬ 
ent interfacial bonding energy with the surrounding. One 
important outcome revealed in Fig. 2(d) is that AEd 
decreases for larger dopants. Therefore, large dopants 
actually have lower site preference, hence are more dif¬ 
ficult to achieve stoichiometric doping. This invalidates 
the previous speculatiorP^. 

One more problem for the large dopants is identified af¬ 
ter fully relaxing the lattice volume and geometry. While 
interkagome doping maintains the original lattice symme¬ 
try, the in-plane doping distorts the lattice from hexag¬ 
onal to triclinic, breaking the perfect kagome plane and 
lowering the in-plane doping energy. Consequently, AEd 
becomes smaller. This effect becomes very significant for 
larger dopants. For example, for CdP^ and Ca^~^, AEd 
decreases from 0.20 and 0.50 eV to -0.12 and 0.11 eV, 
respectively. It means that upon doping, a large frac¬ 
tion of dopants will substitute the in-plane sites, which 
in turn distorts the lattice and hinders stoichiometric se¬ 
lective doping. In contrast, for and which 

have similar radius to Cu^~^, this problem does not occur: 
AEd decrease from 0.76 and 0.50 eV to 0.72 and 0.46 eV, 
respectively, which are still sufficiently large to suppress 
in-plane doping. Therefore, we conclude that for our pur¬ 
pose and Zri^~^ are the most ideal dopants. 

To further examine whether the selective doping of 
Mg/Zn at the interkagome sites can sustain up to 
the stoichiometric limit to form Cu^Mg{OH)QFBr and 
Cu 3 Zn{OH)QFBr compounds, we proceed by progres¬ 
sively increasing the amount of interkagome dopants. Ed 
is found to be nearly independent of the doping concen¬ 
tration (xd)- This property indicates that the interac¬ 
tion between the dopants is weak, which is important for 
reaching the the stoichiometric limit. Otherwise, dopants 
may form clusters or hinder further doping process. 
We have also checked that the lattice and Cu kagome 
planes in Cu 3 Mg{OE[)^FBr and Cu 3 Zn{OE[)QFBr re¬ 
main stable under structural relaxation. 

It is worth making some comparison to Herbert¬ 
smithite within the present methodology. The first crit¬ 
ical issue is the degree of equilibrium anti-site disorder 
in the two systems. The effect of Cu/Zn anti-site disor¬ 
der in Herbertsmithite has been a long-lasting debat^i^. 
Such disorder is inevitable in doped Barlowite as well, 
because it leads to an increase of the configuration en¬ 
tropy. Under constant temperature and volume, the 
equilibrium is reached by minimizing the free energy 
F{xa) = E{xa) — TS{xa) with respect to the anti-site 
concentration Xa- Without loss of generality, it is conve¬ 
nient to set E(0)=S(0)=0 for the stoichiometric systems. 
Accordingly, when Xa pairs of Mg(Zn) and Cu per em- 
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pirical formula switch sites, the energy increase is sim¬ 
ply E{xa) = Xa^Ed, given that the interaction between 
the dopants is weak. The entropy increase per empirical 
formula can be analytically derived as S(xa) = ks InOo, 
where ^la is the anti-site configuration number. With the 
aid of the Sterling’s approximation, the final form of en¬ 
tropy is S{xa) = -fcBln[(3 - Xaf~'^‘^{l - 
The minimal point of F{xa) is calculated numerically, 
which can be expressed as a function of AEd/ksT [Fig. 
2(c)]. By using the calculated values of AEd (Mg: 0.72 
eV; Zn: 0.46 eV) and the experimental growth tempera¬ 
ture (T=393K)I^, the anti-site disorder in Mg/Zn doped 
Barlowite is predicted to be below 0.1%. For compari¬ 
son, we have also calculated AEd = 0.30eF for Herbert- 
smithite and with T=483KI^, the equilibrium disorder 
is calculated to be 5%, comparable to the experimental 
estimatiorP^. Therefore, the degree of anti-site disorder 
in the doped Barlowite is expected to be at least one 
order of magnitude lower than that in Herbertsmithite, 
owing to a higher AEd as well as a lower growth temper¬ 
ature. This difference can be significant to help clarify 
the effects of disorder on the QSL phase. 

Secondly, we do a comparison of single-electron band 
structure. Figure 3 shows the band structures of Hebert- 
smithite, undoped and doped (xd = 1) Barlowite marked 
with orbital compositions. Despite the absence of strong- 
correlation effects, the DFT single-electron band struc¬ 
ture properly describes single-electron hopping processes, 
which serve as the guide to the AFM superexchange. For 
Herbertsmithite, there is a set of bands around the Fermi 
level (between -0.5 eV and 0.75 eV), gapped from the un¬ 
derlying valence band^^ . These bands primarily arise 
from Cu [blue cross in Fig.3(a)] and the adjacent O (not 
shown), exhibiting the typical features of NN hopping on 
a 2D kagome lattic^^. For Barlowite [Fig.3(b)], around 
0.5 eV the band dispersion is similar to that in Herber- 
smithite with the band composition primarily from the 


in-plane Cu, indicating similar hopping amplitude within 
the kagome planes. This is in agreement with the exper¬ 
imental fact that the Curie-Weiss constants for Herber- 
smithite and Barlowite are clos^^. However, around the 
Fermi level, interkagome Cu not only contributes extra 
bands, but also strongly mix with the Cul bands. This 
result suggests considerable coupling b etwee n Cul and 
Cu2, as pointed out by previous studied^oiiu effect 
of replacing Cu2 with Zn or Mg is remarkable [Figs.3(c) 
and (d)]. After doping, the complexities of interkagome 
coupling are removed. Both Zn and Mg states are far 
from the Fermi level, leaving clean Cul bands around 
the Fermi level. The overall band dispersion also becomes 
closer to Herbertsmithite. The energy states below the 
Fermi level contain contribution from the halogen atoms, 
i.e. Cl and Br in Herbertsmithite and Barlowite, respec¬ 
tively. We note that both Cl~ and Br~ are spin zero 
and far away from the superexchange path between Cul 
ions. Therefore, these orbitals do not play an important 
role in the magnetic properties of the materials. 

In conclusion, based on the DFT calculations, we 
identify Cu^Mg{OH)QFBr and Cu'iZn(OH)(,FBr as 
the most promising targets to realize the stoichiomet¬ 
ric doped Barlowite. The distinct advantages include no 
lattice distortion, high site selectivity and low anti-site 
disorder. The standard doping energy is comparable to 
(for Zn) or even lower (for Mg) than that for growing 
Herbertsmithite. Therefore, these targets may be read¬ 
ily synthesized using similar experimental conditions as 
used for Herbertsmithite. The remaining open question 
is how the doped Barlowite behaves magnetically under 
low temperature: will it be tuned into the same phase 
as Herbertsmihite or stay as the undoped Barlowite. For 
either case, the effective doping of this new material as 
we propose here serves as a useful guide to future exper¬ 
iments in a pursuit to reveal key factors towards QSL. 
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